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Positioned between crystalline solids and liquids, disordered many-particle systems which are 
stealthy and hyperuniform represent new states of matter that are endowed with novel physical 
and thermodynamic properties. Such stealthy and hyperuniform states are unique in that they are 
transparent to radiation for a range of wavenumbers around the origin. In this work, we employ 
recently developed inverse statistical-mechanical methods, which seek to obtain the optimal set of 
interactions that will spontaneously produce a targeted structure or configuration as a unique ground 
state, to investigate the spin-spin interaction potentials required to stabilize disordered stealthy 
hyperuniform one-dimensional (ID) Ising-like spin chains. By performing an exhaustive search 
over the spin configurations that can be enumerated on periodic ID integer lattices containing N = 
2, 3,..., 36 sites, we were able to identify and structurally characterize all stealthy hyperuniform spin 
chains in this range of system sizes. Within this pool of stealthy hyperuniform spin configurations, 
we then utilized such inverse optimization techniques to demonstrate that stealthy hyperuniform 
spin chains can be realized as either unique or degenerate disordered ground states of radial long- 
ranged (relative to the spin chain length) spin-spin interactions. Such exotic ground states are 
distinctly different from spin glasses in both their inherent structural properties and the nature of 
the spin-spin interactions required to stabilize them. As such, the implications and significance of 
the existence of such disordered stealthy hyperuniform ground state spin systems warrants further 
study, including whether their bulk physical properties and excited states, like their many-particle 
system counterparts, are singularly remarkable, and can be experimentally realized. 


I. INTRODUCTION 

There has been a tremendous amount of interest re¬ 
cently in the creation of “materials by design,” that is, 
the directed and systematic search for new materials 
that possess prescribed desirable properties. Computa¬ 
tional methods developed for this purpose will play a cen¬ 
tral role toward this goal. Inverse statistical-mechanical 
methods allow for a new mode of thinking about the 
structure and physical properties of condensed phases of 
matter,^ and are ideally suited for materials discovery by 
design. 

Much of statistical mechanics centers around finding 
the structural and bulk physical properties for a given 
many-particle (or many-spin) system Hamiltonian, what 
we refer to as the “forward” problem of statistical me¬ 
chanics.In this regard, so-called inverse statistical- 
mechanical methods have been devised that yield opti¬ 


mized interactions that robustly and spontaneously lead 
to a targeted many-particle configuration of the system or 
targeted set of physical properties for a wide range of con¬ 
ditions. An interesting class of target configurations 
that have been studied are classical many-particle ground 
states of varying complexity and novelty.^’®^^^ One can 
also design the system to have exotic physical proper¬ 
ties, such a negative Poisson ratio and negative thermal 
expansion coefficients over a range of temperatures. 

Recently, we have generalized inverse statistical- 
mechanical methods to the case of two-state Ising spin 
systems with radial spin-spin interactions of finite range 
(i.e., extending beyond nearest-neighbor sites).Our 
interest in these initial studies was to find the optimal 
set of shortest-range pair interactions whose correspond¬ 
ing ground state was a targeted spin configuration. The 
possible outcomes for a given target configuration were 
classified into whether or not the targeted ground-state 
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spin configurations were unique, degenerate, or neither. 
In general, inverse techniques applied to targeted spin 
configurations could have implications for the design of 
solids with novel magnetic and electronic properties, e.g. 
see Refs. 21 and 23. 

Our objective in this work is to determine whether 
one-dimensional (ID) spin chains having target configu¬ 
rations that are disordered, stealthy, and hyperuniform 
can be made to be classical ground states. This choice is 
motivated by the fact that such exotic amorphous states 
of matter in the context of particle systems offer fascinat¬ 
ing open theoretical questions and possess novel physi¬ 
cal properties. Hyperuniform many-particle systems pos¬ 
sess anomalously suppressed infinite-wavelength density 
fluctuations, as quantified by the number variance. 
Disordered hyperuniformity occurs in a variety of physi¬ 
cal systems, including large-scale structure in the Uni¬ 
verse,the arrangement of avian photoreceptors,^® 
driven non-equilibrium systems,^®’®® dynamics in cold 
atoms,surface-enhanced Raman spectroscopy,®^ tera¬ 
hertz quantum cascade lasers,®® wave dynamics in disor¬ 
dered potentials based on super-symmetry,®"*’ and certain 
Coulombic systems.®® 

Stealthy hyperuniform systems are those that com¬ 
pletely suppresses single scattering for a range of long 
wavelengths.®®’®^ It has been shown that systems of par¬ 
ticles interacting with certain long-ranged (ie., on the 
order of the chain length) pair potentials can counterin¬ 
tuitively freeze into classical ground states that are disor¬ 
dered stealthy hyperuniform. By mapping such stealthy 
configurations of particles into network solids, the first 
disordered cellular solids with complete isotropic pho¬ 
tonic band gaps comparable in size to photonic crystals 
were discovered.®®’®® 

Here we begin a program to apply inverse techniques 
to explore whether ID spin chains that are disordered, 
stealthy and hyperuniform can be made to be classi¬ 
cal ground states of radial spin-spin interactions. Such 
spin chains would suppress single scattering for a range 
of wavenumbers around the origin, implying anoma¬ 
lously suppressed magnetization fluctations at long wave¬ 
lengths. By focusing only on configurations that are 
stealthy at the smallest positive wavenumber, we are 
able to enumerate all periodic stealthy spin configura¬ 
tions containing N = 2, 3,..., 36 spins on the ID inte¬ 
ger lattice under periodic boundary conditions and then 
structurally characterize these configurations by comput¬ 
ing pair correlation functions, structure factors, and de¬ 
gree of disorder. We select from them the disordered con¬ 
figurations and determine via inverse optimization tech¬ 
niques whether they can be made to be ground states. We 
discover that stealthy hyperuniform spin configurations 
can be realized as either unique or degenerate disordered 
ground states of radial long-ranged (*.e., relative to the 
chain length) spin-spin interactions. Such exotic ground 
states are distinctly different from spin glasses^®^^^ in 
both their inherent structural properties and the nature 
of the spin-spin interactions required to stabilize them 


(e.g., these disordered stealthy hyperuniform spin sys¬ 
tems exhibit no single-scattering for large wavelengths 
and are stabilized by deterministic rather than stochas¬ 
tic spin-spin interactions). 

Section H briefly reviews basic concepts, including def¬ 
initions of stealthiness, hyperuniformity, and a disorder 
metric for ID spin chains. Section HI describes the 
inverse statistical mechanics procedure used to gener¬ 
ate stabilizing radial spin-spin interactions. Section IV 
presents our major results, including an enumeration 
of stealthy hyperuniform spin configurations containing 
N = 2, 3,..., 36 spins on the ID integer lattice under 
periodic boundary conditions, their characteristics, and 
their ability to be spontaneously generated as ground 
states by the inverse methodology. Finally, in Sec. V, 
we close with concluding remarks. 

II. THEORETICAL BACKGROUND 

We are interested in a generalized class of spin Hamil¬ 
tonians in arbitrary Euclidean space dimensions. In this 
initial study, we focus on the spin-spin interaction Hamil¬ 
tonian of spin chains on the ID integer lattice given by 
the distance-dependent version of the Ising model 

H = -J2 ( 1 ) 

R 

where —1 < J(i?) < 1 can in general be a long- 
ranged spin-spin interaction potential with a radial ex¬ 
tent on the order of the chain length (i? « N). Such 
long-ranged spin-spin interactions have been studied by 
various investigators using standard forward statistical- 
mechanical techniques.The inverse statistical me¬ 
chanics methodology involves optimally tuning the J{R) 
parameters such that the targeted spin configuration is 
spontaneously produced as a unique ground state accord¬ 
ing to the criteria and procedures described in Sec. HI. 
In what follows, we briefly review the basic concepts 
and definitions that will be used throughout the pa¬ 
per, including collective coordinates, hyperuniformity 
and stealthiness, and order metrics. 

A. Collective coordinates 

The configurations we designate as targets in the in¬ 
verse methodology are defined through collective density 
coordinates. Specifically, we focus on ID spin chains 
composed of N spins aj = ±1 positioned on the sites 
of the integer lattice Rj = j in a fundamental unit cell 
under periodic boundary conditions. The spin collective 
density variable of the chains is the Fourier transform of 
the spin density 

N 

p,{k)=Y^a,e^\ (2) 

i=i 
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which can be viewed as a transformation from the fi¬ 
nite set of spin configurations to the complex functions 
Pa{k) that depends on the infinite set of wavenumbers 
k on the integer lattice in reciprocal space.The pos¬ 
sible wavevectors for a chain on the integer lattice are 
k = 2'Kn/N for integer n. The spin structure factor is 
then defined as 


S{k) = ^pa{k)pa{-k), (3) 

which is a real function with inversion symmetry. For 
spins on the integer lattice, the structure factor is peri¬ 
odic with S{k) = S{k + 27r). This symmetry property 
combined with Eq. (3) imply that S{k) is symmetric 
about k = TT for integer spin chains. In addition, through 
Fourier inversion of the structure factor we can obtain 
the radial spin-spin correlation function 
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We note here that the structure factor is expressible in 
terms of the spin-spin correlation function: 

N 

S{k)=Y,S2{r)e^'‘"- 

r—1 

B. Hyperuniformity and stealthiness 

The spin configurations of interest in this study possess 
both hyperuniformity and stealthiness. Hyperuniform 
spin systems possess anomalously suppressed inhnite- 
wavelength magnetization fluctuations, that is^^ 

lim S(k) = 0, 

k^O 

a definition that explicitly omits forward scattering. 

All periodic patterns are hyperuniform, but disordered 
configurations, such as liquids and structural glasses, gen¬ 
erally are not. It is well known that the k —>■ 0 limit of the 
structure factor for particle systems in thermal equilib¬ 
rium is proportional to the isothermal compressibility. 

Stealthy configurations are those in which the struc¬ 
ture factor is zero for a range of wavevectors. In the 
special case in which iS'(fc) = 0 for 0 < fc < AT (where 
AT > 0 is called the exclusion zone radius), the stealthy 
configurations are also hyperuniform, which is the case of 
interest in the present paper. Stealthy patterns are char¬ 
acterized by long-ranged correlations due to the suppres¬ 
sion of large-scale density fluctuations. In this sense, dis¬ 
ordered stealthy hyperuniform patterns have a “hidden 
order” on large length scales.There has been a recent 
effort to generate stealthy particle configurations through 
“stealthy” potentials; see Ref. 37 and references therein. 
These stealthy potentials are radial pairwise interaction 
potentials restricted to a range of wavevectors. These 


FIG. 1. (a) A visualization of how p{k) cancels at the smallest 
positive wavenumber k = jN for (b) the depicted stealthy 
spin configuration of size A = 12 on the ID integer lattice. 
The unit vectors represent the positive spin [oj = -|-1) ex¬ 
ponential terms from the collective density variables 

p(k = 27:jN) = The negative spin terms 

also cancel in a similar fashion. The plot contains a gray 
regular polygon including all possible unit vectors. The unit 
vectors that are summed for the positive spins in this partic¬ 
ular configuration are colored according to how they cancel 
with other vectors. In this case, the scattering cancellation 
can be decomposed into two doublets (green and blue) and 
one triplet (red), which correspond to the indicated spins in 
(b). 


potentials have been shown to produce highly degener¬ 
ate disordered, hyperuniform, stealthy configurations.^^ 
In this study, we explore the discrete spin chain analogs 
of these continuous stealthy hyperuniform particle sys¬ 
tems via the inverse statistical mechanics methodology. 
Specifically, we investigate periodic ID spin chains on the 
integer lattice that are both stealthy and hyperuniform, 
so that S(k) =0 for at least the range of wavenumbers 
0 < k < K = 2 tt/N. Obtaining stealthy spin chains for 
AT > 0 can be non-trivial, as demonstrated in Fig. 1. The 
complex exponentials in the collective density variables 
[cf. Eq. (2)] must cancel at least at the smallest posi¬ 
tive wavenumber k = 2 tt jN for the spin configuration in 
question to be considered stealthy hyperuniform. 


C. Order metric 

We are particularly interested in stealthy hyperuniform 
configurations that possess a high degree of disorder. As 
discussed in the Appendix, periodic spin configurations 
can appear “reducible.” Eor reducible stealthy spin pat¬ 
terns, the unit cell of the pattern can be expressed as a 
repetition of a smaller unit cell. These patterns will of 
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course appear to be more ordered. Of greater interest to 
us are those spin patterns that cannot be decomposed 
into smaller unit cells: the irreducible patterns. These 
patterns are more disordered. 

We quantify the degree of order between configurations 
of the same size N and magnetization (ct) = 
by using the order metric defined by Ref. 37 and used for 
two-dimensional stealthy spin systems in Ref. 51: 

r = ^ [S'(/c) - S'o(/c)]^ (5) 

k 

where the summation is over k = 
27r/A^, 47r/7V,..., 27r(7V — l)/iV and So{k) is a refer¬ 
ence structure factor. For our purposes, we set this to a 
constant value 

So{k) = l-{af, (6) 

which is the expected structure factor of an ensemble 
of uncorrelated Poisson spin patterns of period N and 
average magnetization (a). Each configuration in the 
Poisson ensemble is generated by setting each spin to 
CTi = -|-1 with probability (cr) and to Ci = —1 otherwise. 

The order metric (5) is meaningful for relative com¬ 
parisons between spin configurations of fixed size N and 
fixed magnetization (cr). Four examples of spin config¬ 
urations with (cr) = 0 and different r are displayed in 
Fig. 2(a)-(d) along with their corresponding structure 
factors. Configurations (a)-(c) are stealthy hyperuni¬ 
form, while (d) is a random Poisson pattern. As one 
can see, configuration (a) is periodic on a unit cell of size 
2 and has a single Bragg peak in 0 < fc < 27r. This large 
deviation of S{k) from 5'o = 1 at fc = tt leads to a large r 
value. Configuration (b), the second most ordered con¬ 
figuration, is almost reflectionally symmetric about the 
center of its unit cell, which leads to two peaks smaller in 
size than in (a). This explains the elevated r value. The 
remaining two configurations (c) and (d) do not have any 
easily discernible order and both have S{k) vary about 
the target structure factor = 1. Yet, we see that the 
order metric can pick up the “hidden order” of (c) the 
stealthy hyperuniform configuration compared to (d) the 
Poisson pattern. 

We note in passing that a direct-space representation 
of the order metric defined in (5) is trivially obtained 
by Parseval’s theorem,which is effectively tantamount 
to replacing S'(fc) with the direct-space two-point func¬ 
tion S 2 {r) [defined by (4)] and summing over r. Such 
direct-space sums have been studied in a variety of other 
contexts (not as order metrics), including spin Hamilto¬ 
nians®^ and as optimization objective functionals to re¬ 
construct digitized binary heterogeneous media.®® 

III. INVERSE OPTIMIZATION TECHNIQUES 

We systematically study all stealthy hyperuniform spin 
configurations under periodic boundary conditions with 



FIG. 2. Four periodic spin configurations with N — 30 sites 
and (cr) = 0. The order metric r of each configuration is dis¬ 
played, followed by the spin configuration, and the structure 
factor, (a)-(c) The top three are stealthy hyperuniform, (d) 
The bottom configuration is a randomly generated Poisson 
spin pattern, which is hyperuniform but not stealthy. 


N = 2,..., 36. Specifically, we enumerate the 2^ possi¬ 
ble configurations and keep only those with a vanish¬ 
ing structure factor for at least the smallest positive 
wavenumber fc = 27t/N. In our enumeration, two con¬ 
figurations that can be produced by translations or re¬ 
flections of each other are considered equivalent and are 
not counted as distinct configurations. 

We characterize the structural properties of the discov¬ 
ered stealthy spin systems, including their degree of dis¬ 
order, and attempt to generate spin-spin interaction po¬ 
tentials, which produce the stealthy spin configurations 
as ground states if such potentials exist. We accomplish 
the latter via the inverse statistical mechanics procedure 
for spin systems developed in Ref. 21. The procedure 
searches for interactions J{R) in two-state spin Hamilto- 
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nians of the form 

N Rc 

^ = -EE J{R)aiai+R 

i=l 

Rc 

= -NY,J{R)S2{R) (7) 

R=1 

with interactions up to a cutoff distance Rc between 
spins. The procedure attempts to produce a target con¬ 
figuration T as, a ground state by maximizing the differ¬ 
ence in energy density (e = H/N) between T and the 
closest energetic competitor C: 

Rc 

Ae^ = e^ = J{R) [S^{R) - Sj (i?)] 

R=1 

subject to the constraint that Ae^ > 0 and — 1 < J{R) < 
1. Note that Rc A Rmax must hold, for a certain maxi¬ 
mum distance Rmax, so that the inverse problem is fea¬ 
sible. 

There are three possible outcomes for the inverse 
method, which can be organized into the following so¬ 
lution classes: 

• Class I: Solutions in which a spin-spin potential 
J{R) is found that generates the target configura¬ 
tion as the unique ground state up to translations, 
reflections, and spin inversion operations. 

• Class II: Solutions in which a spin-spin potential 
is found that generates the target configuration as 
a non-unique ground state, degenerate in spin-spin 
correlation 52 (i?) [cf. Eq. (4)] with other spin con¬ 
figurations. 

• Class III: Solutions that are neither class I nor II. 

For a given competitor list, this optimization problem 
is solved using linear programming (to within numerical 
precision) and then simulated annealing with Metropo¬ 
lis Markov Chain Monte Carlo sampling is used to gen¬ 
erate additional competitor spin configurations. Note 
that Class III solutions only arise when the optimization 
is deemed unfeasible by the linear programming meth¬ 
ods. As such, there are no false-positive Class III so¬ 
lutions. For more details on the inverse statistical me¬ 
chanics method applied to spin systems, see Refs. 21 and 
22 . 


IV. RESULTS AND DISCUSSION 

With a complete list of stealthy spin configurations of 
size V = 2,..., 36 at our disposal, we now examine their 
statistical properties, such as magnetization and degree 
of disorder, and classify their “groundstateability,” i.e., 
whether such target spin configurations can be stabilized 
as unique or degenerate ground states. 


A. Characterization of enumerated configurations 

1. Spin configurations of fixed size and magnetization 

The enumerated stealthy spin configurations can be 
expressed in terms of spin configurations of fixed size 
and magnetization, two defining properties of generic 
spin patterns. Of primary interest is how the number 
of stealthy spins grows with system size. Figure 3(a) dis¬ 
plays the abundance of stealthy hyperuniform spin con¬ 
figurations relative to all possible 2^ configurations. The 
general trend for stealthy spin configurations of size N is 
exponential growth in the number of configurations 
at a rate of a » 1.4 < 2. Spin configurations with prime 
N have only the two trivial stealthy configurations of all 
up spins and all down spins. Moreover, most stealthy 
spin configurations of size TV < 36 are reducible to fun¬ 
damental unit cells of size less than N. As discussed in 
the Appendix, this property comes from the fact that all 
stealthy spin configurations oi N = pq^ where p and q are 
prime, are reducible to smaller bases. The values of N 
that do not follow this pattern admit irreducible stealthy 
configurations. They are shown in red in Fig. 3(a) and 
are labeled “irreducible.” 

The probability distribution of magnetization for con¬ 
figurations of stealthy spins of size N, shown in Fig. 3(b), 
is symmetric and binomial-like, though with some key 
differences. For finite system sizes, depending on TV, a 
spin configuration distribution is peaked at either (a) = 0 
or (cr) = ±1/N. This behavior is determined by the 
prime factorization of N and the abundance of reducible 
configurations of size N. For comparison. Fig. 3(b) shows 
the magnetization distribution of all spin configurations 
on chains of length N, which exactly follows a binomial 
distribution. It is interesting to note that the distribu¬ 
tions of magnetization of the stealthy configurations tend 
to oscillate about the binomial profile. This seems to be 
a finite size effect that has to do with the fact that more 
stealthy configurations with an even number of up spins 
exist than with an odd number. 


2. Order and disorder in stealthy spin configurations 

Our enumeration reveals a range of disorder for hype¬ 
runiform stealthy spin systems. The distributions of dis¬ 
order for stealthy spin configurations, according to the 
order metric defined in Eq. (5), are displayed in Fig. 4. 
Figure 4(a) shows the distribution of the order metric r 
among irreducible only and all stealthy spin configura¬ 
tions of size TV = 30 and (cr) = 0. The distribution of t 
tends to spread out as the system size increases and is 
positively skewed. Figure 4(b) shows the distribution of r 
for N = 36 and (cr) = 0 stealthy configurations and high¬ 
lights the importance of irreducibility in stealthy spin 
systems. For size N = 36, the irreducible configuration 
are the most disordered of all the stealthy configurations. 
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FIG. 3. (a) The number of stealthy hyperuniform spin conhgurations for spin configurations of periodic unit cell sizes between 
N = 2 and 36. In gray, for comparison, is the line 2^. In black are all of the stealthy hyperuniform configurations. In red 
are the irreducible stealthy configurations, which cannot be expressed in unit cells smaller than N. (b) The distribution of 
magnetization for irreducible spin configurations. For comparison, the symmetric binomial distributions about (a) — 0 are 
shown as dashed lines. 



0.05 

0.04 

0.03 

0.02 

0.01 

0 


---All 

— Irreducible 


H M 




(a) 


'0 50 100 150 200 

X 

(b) 


FIG. 4. Histograms representing the probability density of the order metric r for stealthy hyperuniform spin chains of size (a) 
N = 30 and (b) N = 3G with zero magnetization. Displayed in dashed red lines are histograms of the r distribution among all 
stealthy patterns of the specified sizes, including the reducible configurations which can be expressed in terms of smaller unit 
cells. A few extremely ordered configurations with r > 500 are not shown. Displayed in solid black lines are histograms of the 
T distribution for irreducible stealthy configurations. Note that there are 2306 size N = 36 stealthy configuration, 872 of which 
are irreducible, and that there are 315 size N = 30 stealthy configurations, 297 of which are irreducible. 


The nature of disorder in our ID spin systems is in¬ 
timately linked to their stealthiness. We measure the 
degree of stealthiness of a material by the number of in¬ 
dependent wavevectors for which S{k) = 0 for k < K, 
which for ID periodic spins chains on the integer lattice 
is M{K) = NK/2t:, proportional to exclusion zone ra¬ 
dius K. The distribution of M{K) for all stealthy spins 
of size = 36 is shown in Fig. 5. We see that the dis¬ 
tribution of M{k) is skewed towards the smaller values 
of K. In general, larger M{K) correspond to configura¬ 


tions with more Bragg peaks in the range of wavenumbers 
0 < /c < 27r and hence larger r values. This is particu¬ 
larly true for the irreducible stealthy spin configurations, 
which only have exclusion radii with M{K) < 2 and are 
the most disordered stealthy configurations. This cor¬ 
relation between stealthiness and order in spin chains, 
along with similar evidence of order in stealthy contin¬ 
uous particle systems with large exclusion zone radii, 
suggests a strong positive relationship between the ab¬ 
sence of magnetic and radiation scattering for large wave- 
































7 



M(K) 


FIG. 5. The distribution of M{K) [number of wavevectors 
for which S{k) = 0 for k < K] for all enumerated stealthy 
configurations of unit cell size N = 36. 

lengths and ordering. Disordered stealthy spin-chains re¬ 
quire a relatively small exclusion zone radius and a nar¬ 
rower range of wavevectors for which scattering does not 
occur. 


B. Groundstateability 

By application of the inverse statistical mechanics 
method developed in Ref. 21, we have been able to show 
that irreducible stealthy hyperuniform spin chains can 
be stabilized as ground states of the long-ranged {i.e., 
relative to the chain length), spin-spin interaction de¬ 
fined in Eq. (7). The “groundstateability” of these irre¬ 
ducible configurations were classified into class I, II, and 
III, according to whether those configurations could be 
uniquely realized as ground states or not, as detailed in 
Sec. III. 

Figure 6 displays the ground-state classification of ir¬ 
reducible stealthy spin configurations. As demonstrated 
in Fig. 4, these configurations are the most disordered of 
the stealthy spin configurations at these sizes. For the 
given stabilizing potential and limited spin chain lengths 
considered here, a majority of stealthy configurations 
are class III. Nonetheless, a surprisingly large fraction 
of stealthy configurations at iV = 18, 24, 30, 36 can be 
uniquely stabilized. Among the stabilized configurations, 
most are class I. Class II stealthy configurations are rare, 
but begin to appear for the larger system sizes. It has 
been demonstrated that, for ID spin chains on the inte¬ 
ger lattice, as N increases the number of S 2 degeneracies 
increases.This trend holds also for stealthy spin chains 
and seems to indicate that more class II configurations 
occur for larger N. This trend of increasing degeneracy 
among stealthy hyperuniform configurations would agree 
with the recent results for the particle analogs. 

We should note that the number of class I and II config¬ 
urations might not be quite accurate due to the heuristic 



(a) 



(b) 

FIG. 6. (a) The number of irreducible stealthy ID spin con¬ 
figurations in classes I, II, and III for sizes N = 12,..., 36 
plotted on a log-scale. The first size with class II configura¬ 
tions is A = 30. A = 12 has one class III configuration. Gon- 
figurations equivalent by a spin-inversion are not counted, (b) 
The relative amounts of these configurations in classes I, II, 
and III. Most of the sizes studied have no irreducible stealthy 
configurations and so are not depicted. Only 10 out of 1068 
irreducible stealthy configurations with period A = 30 and 
52 out of 5362 irreducible stealthy configurations with period 
A = 36 are class II. 


nature of the zero-temperature competitor-based method 
implemented.^^ Class III configurations cannot be incor¬ 
rectly classified, as discussed in Sec. III. However, true 
class III configurations can be incorrectly classified as 
class I or II if the competitor space is not adequately 
sampled by the simulated annealing procedure. Incor¬ 
rect classification is more likely for large fundamental 
unit cells N and large potential potential cutoffs Rc, 
where the solution space becomes large and full of local 
minima. Nonetheless, we took care to verify the class I 
and II configurations found through the inverse method. 
After obtaining a solution, we ran 10 - 50 of iterations of 
simulated annealing optimizations to confirm the target 
configuration as the ground state. Each simulated an¬ 
nealing run started with different initial conditions and 
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FIG. 7. The probability distributions of relative potential 
cutoffs Rc/N, where Rc is an integer between 1 and N, for 
= 18 and N = 36 periodic ID irreducible stealthy spins 
systems on the integer lattice. Larger systems tend to have 
a greater spread in relative potential cutoff centered around 
lower values. For = 18 and A’ = 24 (not shown), all 
potential cutoffs were at least half of the period length Rc > 
N/2. While N = 30 (not shown) and A = 36 contained many 
configurations which could be represented as ground states of 
potentials of length Rc < N/2. 


involved 10,000 - 30,000 Monte Carlo sweeps. 

The J(i?) potentials engineered by the inverse method 
are the shortest possible stabilizing interactions. Inter¬ 
estingly, these interactions are long-ranged, relative to 
the chain length A, and display a surprising range of 
cutoff distances, with 9 < Rc < 23 for size A = 36 
stealthy spin chains. Figure 7 shows the distributions of 
potential cutoffs Rc for Class I and Class II irreducible 
stealthy spin configurations of size A = 18 and A = 36. 
An important characteristic of these distributions is that 
the variance in the range of the potential increases with 
A. This is promising for the future engineering of disor¬ 
dered stealthy spin configurations. It would seem plausi¬ 
ble for short-range {Rc/N 1) interactions to produce 
disorder, stealthiness, and hyperuniformity. 

There is a large variety of interactions that can stabi¬ 
lize stealthy hyperuniform spin chains as ground states. 
To showcase this variety in size A = 36 spin chains, 
we pick three potentials with Rc from different points 
{Rc = 9,12,14) in the distribution depicted in Fig. 7. 
Figure 8 shows the three stabilizing potentials, along 
with the class I stealthy spin configurations they gen¬ 
erate and the spin-spin correlation of the configurations. 
Figure 8(a) shows a (cr) = 0 conhguration generated by 
a short {Rc = 9) potential. Both the spin-spin correla¬ 
tion and the potential J{R) vary greatly on short length 
scales. Figure 8(b) shows a longer {Rc = 12) potential, 
which, along with the spin-spin correlation, changes more 
smoothly as a function of distance. The longest potential 
displayed {Rc = 14) in Figure 8(c) shows the smoothest 



FIG. 8. Three representative stealthy class I spin configura¬ 
tions of size 30 with (a) = 0. Depicted are their unique spin- 
spin correlations S 2 {R) in blue and their stabilizing potentials 
J{R) in red with potential cutoffs {a.)Rc = 9, (b) Rc = 12, 
and (c) Rc = 14. Note that the J{R) are not unique and that 
other stabilizing interactions could exist, though not with cut¬ 
offs less than the given de¬ 


stabilizing potential and spin-correlation of the three con¬ 
figurations shown. These three class I spin configurations 
and their corresponding stabilizing potentials are repre¬ 
sentative. 

Our characterization of class II configurations is more 
limited. We have only managed to observe a few dozen 
such configurations for size 30 and 36 spin chains. All size 
30 class II stealthy configurations are shown in Fig. 9. 
For size 30 and 36 stealthy hyperuniform spin chains, 
all class II configurations came in S' 2 -degenerate pairs, 
rather than in larger degenerate groups. It seems likely 
that for larger A, degeneracy should increase given the 
large combinatorial increase in the number of possible 
configurations and ways to construct stealthy, hyperuni¬ 
form patterns. 
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FIG. 9. All class II stealthy spin configurations of size 30. Five pairs of S 2 (i?)-degenerate configurations are grouped horizontally. 


V. CONCLUSIONS 

In this work, we have shown that disordered stealthy 
hyperuniform spin configurations can be realized as ei¬ 
ther unique or degenerate ground states of radial long- 
ranged (relative to the chain length) spin-spin interac¬ 
tions. Stealthy hyperuniform states are unique in that 
they are transparent to radiation for a range of wavenum¬ 
bers around the origin, implying anomalously suppressed 
magnetization fluctuations at long wavelengths. The 
discovered exotic disordered spin ground states, dis¬ 
tinctly different from spin glasses, are the spin analogs of 
disordered stealthy hyperuniform many-particle ground 
states^®’^^ that have been shown to be endowed with 
novel photonic properties.Thus, stealthy hyper¬ 
uniform spin systems offer potentially exciting new av¬ 
enues for future research, as we will elaborate below. 

It is useful to summarize how we came to ascertain 
that such disordered spin ground states exist. First, 
we performed an exhaustive enumeration over the en¬ 
tire set of spin configurations that exist on periodic ID 
integer lattices containing N = 2, 3,..., 36 sites in or¬ 
der to identify and characterize all stealthy hyperuni¬ 
form spin chains in this system size range. In doing so, 
we found that the number of stealthy hyperuniform spin 
configurations grows exponentially with system size (i.e., 
~ 1.4^), implying that the fraction of these unique spin 
configurations goes to zero in the thermodynamic limit. 
Furthermore, the distributions of magnetization and or¬ 
der in these stealthy hyperuniform spin configurations 
show strong deviations from the binomial distributions 
of these quantities that are characteristic of the set of all 
enumerated spin configurations. 

To study the “groundstateability” of disordered 
stealthy hyperuniform spin configurations, we employed 
recently developed inverse statistical mechanics tech- 
niques^^’^^ in conjunction with a class of Hamiltonians 
that allow for radial pairwise spin-spin interactions that 
extend well beyond nearest-neighbor lattice sites. Al¬ 
though many of these spin configurations cannot be sta¬ 
bilized as unique (Class I) or S' 2 (i?)-degenerate (Class II) 
solutions within this set of allowed spin-spin interaction 
types, we did identify a significant number of Class I and 


II disordered stealthy hyperuniform spin chains, in par¬ 
ticular for the largest lattices considered herein. Interest¬ 
ingly, the spin-spin interaction potentials that were able 
to spontaneously produce these disordered spin chains 
display a wide radial extent, spanning from « 30% to « 
60% of the length of the entire underlying integer lattice. 
Although these long-range interactions (with respect to 
the size of the lattice) tend to be attractive over their 
defined radial extent, there still exists a great variety in 
the shapes and relative magnitudes of these spin-spin in¬ 
teraction potentials. As such, these findings are promis¬ 
ing indicators that such interactions—and therefore such 
exotic spin configurations—can in fact be realized exper¬ 
imentally. 

Our fundamental understanding of disordered stealthy 
hyperuniform spin configurations and the spin-spin inter¬ 
action potentials required to stabilize them is currently 
in its infancy. Hence the implications of the existence of 
disordered stealthy hyperuniform spin ground states es¬ 
tablished herein provide fertile ground for future research 
directions. We are interested in determining whether the 
excited states and bulk physical properties of these sys¬ 
tems are singularly remarkable and characterizing these 
potential novel states of matter. The ground-state classi¬ 
fication of these systems in the thermodynamic limit still 
remains an outstanding open question as well as the ef¬ 
fects of thermal and/or quantum mechanical fluctuations 
on their stealthy and hyperuniformity properties. Such 
knowledge would be invaluable in realizing the rational 
design of these exotic spin systems and exploring their 
potential technological applications. 
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Appendix: Reducibility of stealthy hyperuniform 
spin chains 

Most stealthy hyperuniform spin chains on the integer 
lattice are represented in a larger than necessary funda¬ 
mental unit cell. These so called reducible configurations 
can be defined with a smaller fundamental unit cell, un¬ 
der which they are not stealthy. Consider, for instance, 
the antiferromagnetic configuration depicted at the top 
of Fig. 2. This configuration is not stealthy, i.e. does 
not have S{k) = 0 at the smallest positive wavenumber 
k = 27r/iV, when represented in a periodic fundamental 
unit cell of size N = 2. 

The unusual trend in the number of irreducible 
stealthy configurations as a function of fundamental unit 
cell size N, shown in Fig. 3(a), suggests that the re¬ 
ducibility of stealthy configurations is related to the 
prime factorization of N. While we could not easily iden¬ 
tify the general trend, we were able to characterize the 
reducibility of a particularly simple factorization of N. 

Conjecture: A spin chain that is stealthy hyperuniform 
on the integer lattice with fundamental unit cell size N = 
pq, where p and q are prime, is always reducible, i.e. can 
be equivalently represented by a smaller unit cell of size 
TV' < TV in which it does not appear stealthy. 

Discussion: We use the notation a — (cti, (T 2 , ..., ctat) 
to signify spin configurations. At size TV, one can eas¬ 
ily construct two crystalline configurations and 
that cause the collective density variables, and hence the 
structure factor, to vanish when k = 2tt/N. First, we 
will show this for the configuration 




(p) 


j = p, 2p,... 
otherwise 


qp 


(A.l) 


Its collective density variables at the smallest positive 
wavenumber can be expressed as 


N 


^(p) 


(2^/TV) = ^ 

i — 1 

^(P) gi{2irj/N) 

(A.2) 

q 

N 


= E 

gi(27mp/Af) _ ^i{2'n'j/N) 

(A.3) 

n—1 

j^np 


II 

(27r/TV) - p^^(27r/TV) 

(A.4) 


We see that the second term can be expressed in terms 
of the first term 


where we used the periodicity of the chain. 

Now we can see that if vanishes, so does As 

(p) 

it turns out, the pj. term does indeed evaluate to zero: 


P^:\27t/N) = 


q 

^i{2-KnplN) 

n—1 


1 


1 — g-i(27rp/W) 

1 

1 — g-i(27rp/W) 
= 0 . 


_ gi(27rp<j/JV)^ 
I'l _ 


(A.7) 

(A.8) 

(A.9) 

(A.IO) 


Therefore, p^P\2Tr/N) = 0 and is stealthy. The same 


logic applies to defined as 


= 


j = 9, 2(7,. 
otherwise 




(A.ll) 


so we can conclude that it is also stealthy. 

Both and are reducible. One might think 
that a new irreducible stealthy configuration a’ could be 
generated by superposing and taking cr' = -|-1 


when = -|-1 or cr)'" = -|-1. Based on the same reason 
presented in Sec. VI of Ref. 37, the superposed config¬ 
uration is stealthy if the sequences {p, 2p ,..., qp} and 
{< 7 , 2q,... ,pg} do not overlap. However, they do overlap 
exactly once at pq. Since the collective density variables 
are linear in their terms, this leaves p'{2tt/N) with one 
term in its summation that does not cancel. Therefore, 
p'{2tt/N) 7 ^ 0 and the overlapped configuration is not 
stealthy. 


Aq) 


P- (27r/TV) = 


,i(2ir(np-|-l)/Ar) 


-k ••• -f 


P-1 

E 

m—l 


g*(2Wiv)p(p)(27r/A) 


Moreover, our enumeration results suggest that a^P'> 
and are the only stealthy configurations that exist for 
N = pq (modulo spin inversion and translation). Both 
are trivial configurations that are reducible to fundamen¬ 
tal unit cells of size N' = p and q respectively. In the 
q reduced TV' cell, these chains have only a single up spin 

gi(27r(np-i-(p-i))/Ar)5y}^jg]^ leads to a term that does not cancel in the density 
n=i p{2'k/N'). Therefore, the reduced configurations are not 

(A.5) stealthy. 

(A.6) Therefore, all stealthy hyperuniform spin chains of size 

N = pq are reducible. 
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